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ABSTRACT 

Context. Dusty debris discs around main sequence stars are thought to be the result of continuous collisional grinding of planetesimals 
in the system. The majority of these systems are unresolved and analysis of the dust properties is limited by the lack of information 
regarding the dust location. 

Aims. The Herschel DUNES key program is observing 133 nearby. Sun-like stars (< 20 pc, FGK spectral type) in a volume limited 
survey to constrain the absolute incidence of cold dust around these stars by detection of far infrared excess emission at flux levels 
comparable to the Edgeworth-Kuiper belt (EKB). 

Methods. We have observed the Sun-like star HD 207129 with Herschel PACS and SPIRE. In all three PACS bands we resolve a ring- 
like structure consistent with scattered light observations. Using a Bootis as a reference point spread function (PSF), we deconvolved 
the images, clearly resolving the inner gap in the disc at both 70 and 100 fim. 

Results. We have resolved the dust-producing planetesimal belt of a debris disc at 100 iJ.m for the first time. We measure the radial 
profile and fractional luminosity of the disc, and compare the values to those of discs around stars of similar age and/or spectral type, 
placing this disc in context of other resolved discs observed by HerschelfDV^ES. 

Key words, stars: individual: HD 207129, stars: circumstellar matter, infrared: stars 



1. Introduction 

Debris discs are composed of dust grains continuously produced 
by the collisional grinding of larger unseen planetesimals . This 
is inferred from the sh ort lifetime of the dust gr ains compared to 
the age of the star (IBackman & Pares"cell 19931) and suggest that 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 
** FNRS Postdoctoral Researcher 



the star around which they are observed has undergone a plan- 
etesimal formation process. In almost all cases, dust production 
is consistent with s teady state attrition of the dust parent bodies 
dLohne et al.l l2008). The vast majority of debris discs are unre- 
solved and we can only use modelling of the disc spectral en- 
ergy distribution (SED) to determine the spatial location of the 
dust and the physical properties of the constituent dust grains. 
In the absence of independent information constraining the spa- 
tial location of the dust, which introduces degeneracies into the 
fitting of the dust grain properties, e.g. between grain size and 
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radial location, a standard dust composition is assumed. In the 
few examples of resolved debris discs, structures that imply the 
presence of a planetary mass body to maintain and sculpt the 
disc are frequently observed, e.g. warps, asymmetries and blob s 
(iGohmowski et al.ll2006l iKalas et alJl206a ISheret et aT]l200l . 
Planets may also be responsible for the inner cavities seen in 
resolved discs and inferred from the SEDs of unresolved ones 
(Wyatt 2008). The resolved disc structure can therefore be used 
as an indirect probe for exoplanets around such stars in regions 
of orbital radius/planetary mass parameter space that are other- 
wise inaccessible to traditional search methods (e.g. radial ve- 
locity, transits) and remain a challenge to direct imaging tech- 
niques. 

HD 207129 ( HIP 107649) was identified as having a de- 
bris disc by IRAS dWalker & Wolstencroft 1988) and followed 
up by both ISO a nd Spitzer ( Jourdain de Muizon et al. 1999t 
iTrilhng et al.ll2008h . Extended emission from the disc has been 
seen in both scattered light (HST) and in thermal infrared 
emission (Spitzer MIPS70, Krist et al. 2dTol). The disc excess 
was also d etected at 160 fim by both iTanner et al.l (l2009h and 
iKristetaU jIolO), though with very diiferent values (155 mJy cf 
250 mJy), illustrating the need for Herschel PACS observations 
to constrain the disc SED. Additionally, the system has been 
identified as a promising candidate for exoplanet searches, due to 
its proximity (1 6 pc), age (1-3 Gyr) and the presence of an inner 
gap in the disc (iJourdain de Muizon et al.l[T999t iBeichman et al.l 
l2010h . 

In this paper we pr esent Herschel (Pilbratt et al. 2010) PACS 
dPoglitsch et al.ll201(3h and SPIRE (.Griffin et al..,201 0) observa- 
tions of HD 207129. The large aperture Herschel telescope pro- 
vides arcsecond resolution allowing detailed imaging of the de- 
bris disc in this system. In addition, greater precision and denser 
coverage of the disc SED can be obtained. Altogether, this al- 
lows better constraints to be placed on the disc's SED and phys- 
ical extent compared to previous observations, thereby allowing 
additional refinement of our models of this solar system ana- 
logue. 



2. Observations and data reduction 

HD 207129 w as obseiTed as pa rt of the DUNES (DUst around 
NEarby Stars; lEiroa et aL I I2OIOI Rodmann et al., in prep.) vol- 
ume limited survey of nearby (d < 20 pc) Sun-like (FGK) stars. 
PACS scan map observations of the star were taken with both 
70/160 and 100/160 channel combinations. Each scan map con- 
sisted of 10 legs of 3' length, with a 4" separation between legs, 
at the medium slew speed (20" per second). The target was ob- 
served at two position angles (70 and 1 10°) in both wavelength 
combinations. SPIRE small map mode observations were also 
carried out on a separate Observation Day (OD) covering a re- 
gion ~ 4' around HD 207129 at the nominal slew speed (30"/s). 
Using five repetitions of the scan map observations reduced the 
expected noise level to close to that of the expected extragalactic 
contribution (7-9 mJy at 250-500 yum) and increased the cover- 
age of the central region of the map allowing pixel sizes smaller 
than the standard values to be used in the image reconstruction 
process (i.e. image scales of 4", 6" and 8" per pixel at 250, 350 
and 500 /urn, limited by the appearance of gaps in the image cov- 
erage near the centre of the map), which was useful for looking 
at extended structure in the source brightness profile. A summary 
of the observations is presented in Table [1] 



Table 1. Summary of Herschel observations of HD 207129. 



Instrument 


Observation ID 


OD 


Wavelengths 


Duration 








Ipm] 


[s] 


PACS 


1342193163/64 


322 


70/160 


276.0 


PACS 


1342193165/66 


322 


100/160 


2250.0 


SPIRE 


1342209300 


544 


250/350/500 


721.0 



Table 2. Summary of Herschel and ancillary infrared/sub-mm 
photometry used in disc SED fitting. Uncertainties l-cr and 
driven by the sky noise. Herschel calibration uncertainties are 
10% for PACS and 15% for SPIRE. 



Wavelength Fo;,, Instrument 

iMin] [mJy] 

9 
18 
24 
32 
70 
160 
160 
70 
100 
160 
250 
350 
500 
870 



2. 1 . Herschel photometry 

PACS data reduction was carried out in HIPE 4.2 (the latest 
available public releaseQ), starting from the level products us- 
ing the standard reduction script. The separate scans at the two 
position angles of each channel pair were mosaiced to produce 
a final image at each wavelength. Image scales for the final mo- 
saics were 1" per pixel for the blue (70/100 yum) images and 
2" per pixel for the red (160 /vm) image. A high-pass filter was 
used to remove large scale background emission from the im- 
ages, with filter widths of 15 and 25" in the blue and red chan- 
nels, respectively. A central region of 30" radius in the images 
was masked from the high pass filter process to prevent the re- 
moval of any faint extended structure near to the source. SPIRE 
data reduction was also carried out in HIPE 4.2, again starting 
from the level data using the standard script and processing 
options. The SPIRE maps were created using the naive scan 
mapper alogrithm. The pixel scales for the SPIRE 250, 350 and 
500 yum images from which photometry was taken were 6", 10" 
and 14". 

PACS fluxes were measured using aperture photometry car- 
ried out using the IDL APER routine. The aperture radius and 
sky annulus dimensions were 20" and 30-40", respectively. 
SPIRE fluxes were measured using an aperture of radius 30" 
(due to the presence of several nearby sub-mm bright back- 
ground objec t s), wh ilst the sky noise values were taken from 
iNguven etal1(l2010l) . 

2.2. Stellar parameters 

HD 207129 (HIP 107649) is a nearby (d ^ 16 + 0.2 pc, 
Ivan LeeuwenI 12007) star with a reported spectral type of G2V 
(from the Hipparcos catalogue. iPerrvman et al.l Fl 9971) or GOV 

' see: http : / /herschel.esac.esa.int/HIPE.download.shtml 



1237 


+ 17 


AKARiaRC PSC, Ishihara et al. (2010) 


263 


± 31 


AKARI/IRC PSC, Ishihara et al. (2010) 


155 


± 5.3 


Spitzer/MIFS, Trilling et al. (2(m). 


111 


± 5.1 


Spitzer/IRS 


278 


± 11 


Spitzer/MlPS, Trilling et al. (2008) 


250 


± 40 


Spitzer/MIPS, Krist et al. (2010) 


158 


± 20 


Spitzer/MIFS. Tanner et al. (2009) 


284 


± 1.5 


Herschel/PACS 


311 


± 1.1 


Herschel/PACS 


211 


± 1.5 


Herschel/PACS 


113 


± 18 


Herschel/SFIRE 


44.3 


± 9 


Herschel/SRIRE 


25.9 


± 8 


Herschel/SRIRE 


5 + 


3 


APEX/LABOCA. Nilsson et al. t2010) 
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Table 3. Physical properties of HD 207129. 



3. Results 



Parameter 


Value 


Distance 


16 + 0.2 pc 


Spectral type and luminosity class 


G2V, GOV 


V,B-V 


5.57, 0.60 mag 


Absolute magnitude My, Bolometric Correction 


4.55, -0.06 


Bolometric luminosity, L, 


1.258 Lo 


Effective temperature 


5912 K 


Surface gravity, log g 


4.44 


Radius, R, 


1.07 Rq 


Metallicity, [Fe/H] 


-0.01 


Rotational velocity, v sin i 


3.71 kms-' 


Rotation period, P 


~ 12.6 days 


Activity, log^Jjj, 


-4.80, -5.02 


X-Ray luminosity, log Lx/L, 


-5.63 


Mass, M, 


1.0-1.15 Me 


Age 


1.5-3.2 Gyr 



(iGrav et al.l l2006h . The bolometric luminosity has been esti- 
mated from the absolut e magnitude a nd bolometric correction 
using measurements by iFlowerl (11996 ). Our adopted values for 
the effective temperature, gravity and metallicity are derived 
from t h e mean of spectroscopic m easur ements fromlSantos et al] 
(12004 . IValenti & Fisciie^ (l2005l) and ISousa et al] (|2008 |). Our 
own estimate of the rotational velocity, from high resolution 
spectra, is u sin / - 3.71 + 1.81 kms"' (Maldo nado e t al., 
in prep.), co nsistent with estimates of .Groot et al.i (1 1996b and 
iTorres et"aD (12006). Assuming the stellar inchnation is that of the 
disc (/ - 60 + 5°, see Sect. 3.2) and the stellar radius is 1.07 R0 
(estimated from the bolometric luminosity and effective temper- 
ature), the rotational period of the star is P = 12.6 days. The star 
is non-a ctive, with measure ments of th e activity ind ex, log R[i^, 
of - 4.8 lHenrv et alld 19961) and - 5.02 lGrav etall (|2006), and a 
ROSAT X-ray luminosity log LxlL^oX - - 5.63. 

The stellar mass estimated from the r adius and gravity i s 
1.15 Mq. From Padova evolutionary tracks (iGirardi et al.ll2002l . 
a mass of 1.0-1.1 Mq and age of ~ 3.2 Gyr a re obtained. 
This age is consistent with the activ ity index of Henry et al. 
( 19961) . using the calibration from iMamaiek & Hillenbrand 
(l2008h . Conversely, the age derived from the rotational period 
using the same calibration is ~ 1.6 Gyr, which is consistent 
wi th the X-ra y luminosity age of 1.8 Gyr using the relationship 
of iGarces et al. (2010), or the results of Giardino et al. (200^ 
(their Fig. 8). Ot her age estimates range from 0.6 (Song et aT] 
120031) to 6.0 Gyr (jLachaume et al.ll999l:IValenti & Fischer.2005F 
iHolmberg etal] 12009). We measure a Lil 6708 A equivalent 
width of 35.5 + 3.3 mA (Maldonado et al., in prep.); this value, 
in conjunction with the estimate derived from the rotation pe- 
riod, points to an age greater than 600 Myr, implying the stellar 
age lies between these extremes. We therefore adopt an age in 
the rang e 1.5-3.2 Gyr as appro priate for HD 207129 (see Figs 
6 and 9. lMaldonado et ani2010l) . The stellar parameters and ob- 
servational properties are summarised in Table |3] 

The stellar photosphere contribution to the total flux was 
computed, using the stellar parameters, from a synthetic stel- 
lar at mosphere model interpol ated from the PHOENIX/GAIA 
grid (iBrott & Hauschild3l2005h . Optical and near infrared pho- 
tometry including Stromgren uvby, Tycho-2 BV and JHK from 
[Aumann & Probst ( 1991) constrain the stellar component of the 
SED, which has b een scaled to the Spit zer IRS spectrum follow- 
ing the method of iBertone et all (l2004 . 



3. 1 . Disc spectral energy distribution 

For the purposes of disc SED modelling, the PACS photome- 
try was supplemented by a broad range of infrared and sub- 
millimetr e observations inclu ding AKARI IRC all-sky survey 
9/18 yum (llshihara et al.ll20 Oft. Spitz.er IRS spe ctrum and MIPS 



24 and 70 A t m pho tometry dTrilling et ani2008l) . S pitzer 160 urn 
dKrist et al.' '201(f), APEX/LABOCA 870 jim (iNilsson et al.1 
|010) fluxes and upper limits from 450/ 850 um JCMT/SCU BA 
(fSheret et al.l 12004)) and 1.2 mm SEST (ISchtitz et al.ll2005h ob- 
servations. The Herschel photometry and all complementary 
data are summarised in Table |2] (the upper limits are not quoted 
in the table because they were not used in the SED fitting pro- 
cess). 

The disc SED, with the Herschel fluxes and ancillary pho- 
tometry can be seen in Fig [1] We have fitted the excess emis- 
sion from HD 207129 with a standard black body model , mod- 
ified beyond 210 jum by a factor of y6 = 1 (see Eqn 6, IWvattI 
2008), from which estimates of the disc fractional luminosity 
and disc orbital radius were derived. The model was fitted to 
the infrared (both Herschel and ancillary) photometry through 
least squares minimisation weighted by the observational uncer- 
tainties. The best fit temperature Tdisc = 50 K and fractional 
luminosity Lir/LboI = 8.3 x 10"^ were obtained with a reduced 
of 1.34. Using the derived stellar parameters and the typical 
disc temperature, we calculate a dust orbital radius of 34 AU 
for black body emission, clearly much smaller than the resolved 
disc which would imply that a disc orbital radius larger than that 
derived from the black body modelling is required. 



100.000 



1.000 



0.100 



0.010 




1.0 10.0 

Wavelength [microns] 



1000.0 



Fig. 1. HD 207129 spectral energy distribution. Red data points 
are the PACS and SPIRE photometry, black data points are 
the ancillary data, including the Spitzer IRS spectrum. Error 
bars are 1 cr and include both calibration uncertainty (10% for 
70/100 yum, 20% for 160 //m) and sky noise. Some error bars 
are smaller than the data point symbol. For wavelengths longer 
than ~ 30 jum the solid line represents a disc model with a 
/lo = 210 /im and P - \, the dashed line represents a disc model 
with a /lo - 63 //m and j3 - 2.7 and the dotted line represent the 
stellar photosphere. 
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Long axis 



Short axis 




5 10 15 

r [arcsec] 



20 



5 10 15 
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Fig. 2. Radial profiles of HD 207129 along the semi-major (left) 
and semi-minor (right) axes of the disc. The red dotted lines are 
the measured profile either side of the disc centre along each 
axis; the solid red error bars mark the mean position at each point 
and the associated uncertainty. The grey line is the deconvolved 
disc profile and the black line is the subtracted stellar profile. 



3.2. Disc images 

The disc of HD 207129 is clearly resolved in the PACS mosaics 
at all three wavelengths and is extended in the two shorter SPIRE 
wavelength maps. The optical position of HD 207129 is (in 
the epoch of the Herschel images) 2l''48"'16.05 ' -47°18'14.47" 
(using proper mot ions from the re-reduction of Hipparcos, 
Ivan Leeuwenll2007l) . In the PACS 70 /im image there is a peak 
in the disc brightness ~ 2" NE from this position, consistent 
with the optical position within the Herschel pointing uncer- 
tainty, with a flux value consistent with the predicted stellar pho- 
tosphere contribution. If this peak were the star, that would im- 
ply that the disc was asymmetric. In the deconvolution presented 
here, we have assumed that the star is at the centre of the disc, 
consistent with the symmetric sc attered light disc observed in 
the HST results (" Krist et al]|2010l) . and that the peak observed in 
the 70 yum image is a dust blob. 

The disc size and position angle were measured by fitting a 
rotated ellipse to the source brightness profile contour of 3 times 
the sky noise value in the PACS images, whilst the disc inclina- 
tion is measured from the ratio of the semi-minor to semi-major 
axes. A summary of these properties is in Tabled The disc ex- 
tent is the same in both the 70 and 100 //m images, though we 
would naively expect that the apparent disc size would increase 



Table 4. Measurement of the disc extent before (left - the 
HWHM of a Gaussian fitted to the source brightness profile) 
and after (right - the radial extent of annulus peak brightness) 
deconvolution, and the inclination and position angle of the disc 
at all three wavelengths. 









Wavelength 






70 yum 


100 yum 


160 yum 


Semi-major axis ["] 


14.5 


9.0 


14.0 


9.0 


17.0 8.5 


Semi-major axis [AU] 


232 


144 


224 


144 


272 136 


Semi-minor axis ["] 


8.5 


4.0 


8.5 


4.5 


12.0 4.0 


Semi-minor axis [AU] 


136 


64 


136 


72 


192.0 64 


Position angle" [°] 


120 




122 




120 


Inclination* [°] 


54 


64 


53 


60 


45 62 



"The disc position angles are the same in both original and 

deconvolved images. 
*The disc minor axis is not resolved at 160 //m in the original images. 



due to the larger beam size and greater contribution from colder 
dust emission. We attribute the similar disc sizes to a combina- 
tion of the similar beam sizes at the two wavelengths (5.7" and 
6.7" FWHM at 70 and 100 yum, respectively) which are both 
small compared to the extent of the disc. In both the 70 //m 
and 100 jum image, there are two peaks in the source brightness 
profile, with the NW peak ~ 10% brighter than the SE one (at 
100 //m). This structure is consistent with the underlying disc be- 
ing a symmetric ring-like structure. Additionally at 70 pm, there 
is a third peak close to the centre of the disc, which is NE (~ 2") 
of the disc centre and stellar position, which may be the result of 
inhomogenities in the disc structure. The disc source brightness 
profile at 160 pm (and in the SPIRE images) is smooth, showing 
no internal structure. 

In order to measure the true extent of the disc, the images 
were deconvolved from the instrument PSF, a technique which 
has previously proved successful (Liseau et al. 2010). An obser- 
vation of a Bootis was used as the PSF model for the deconvo- 
lution, rotated to match the roll angle of the telescope at the time 
of observing HD 207129. 

Deconvolution was a two step process; the stellar photo- 
sphere contribution was removed from each image by subtrac- 
tion of a PSF with a peak scaled to the predicted photospheric 
flux level in that image and centred on the stellar position deter- 
mined from measurement of the 70 pm isophotes. The optical 
and isophote derived positions are in good agreement and the 
small offset between them does not impact upon the findings of 
this analysis. 

After star subtraction, the image was deconvolved using 
three separate methods (modified Wiener, Richardson-Lucy and 
van Cittert) to check the suitability of the individual methods and 
the repeatability of any structure observed in the deconvolved 
images. It was found that all three methods clearly produced a 
ring structure in the deconvolved 70 and 100 yum images, though 
with varying noise patterns. Deconvolution of the 160 yum image 
via modified Wiener and van Cittert methods resulted in images 
with the disc structure as a pair of blobs either side of the stellar 
position. Using the Richardson-Lucy method the disc structure 
at 160 yum was recovered as a broken ring surrounding the stellar 
position with a clear central gap. We do not interpret the clumpy 
blobs in the deconvolved disc images as representing real struc- 
ture in the debris disc. The result of the image deconvolution 
using the modified Wiener method can be seen in Fig. [3] 

The deconvolved radial extent of the disc was measured from 
the position of the peak brightness of the observed disc annulus 
along its major and minor axes. In the deconvolved 100 yum im- 
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age, the measured radial extent was 9" + 2" (144 + 32 AU, see 
Fig|2]l, with similar values from the 70 and 160 /im images. We 
measure the position angle of the disc to be 120 + 5° through 
fitting via a least squares minimisation algorithm. From the ra- 
tio of the semi-major and semi-minor axes, the inclination to be 
54 + 5° (60 ±5° in the deconvolved images) . These results ar e 
consistent with previous HST measurements ( Krist et a l] l20T(l . 
We see no shift in the positional angle of the disc toward longer 
wavelengths. 

3.3. Discussion 

We have directly resolved the structure of a ring-like debris disc 
for the first time in the far infrared. The dis c radial extent of 
~ 140 AU is comparable to that of Fomalhaut dKalas et alj2005h 
or HD 107146 (Ardila et al . 2004). 

The disc radius derived from the standard black body model, 
modified beyond 210 //m by a factor of/3 - 1, was 34 AU. This is 
a large underestimate of the true extent, in direct conflict with the 
resolved images. It serves to illustrate the dangers of modelling 
debris discs using only the SED and/or black body thermal emis- 
sion models, which neither assum e, nor tell yo u anything about, 
the dust grain optical properties (IWvattll2008l) . Using a disc of 
dust grains that emit as modified black bodies and 10 //m dust 
grains, with a break shortward of the SED peak, the disc radius 
derived from the SED matches the observed images but fails to 
reproduce the sub-millimetre slope. The observed dust optical 
and thermal emission properties cannot be reconciled via Mie 
theory, either A disc composed of small (< 1 fim) dust grains 
would, according to Mie theory, reproduce the symmetric scat- 
tered light image, t hough the dust w ould then be too warm to 
reproduce the SED (iKrist et al.ll2010h . 

We measure a dust fractional lumino sity of 8.3x10"^, aro und 



Habing et al.l l2001h and 



six times greater than Vega (1 .5 x 10 

half that of q' Eri (1.5 x lO""*, iLiseau et al.ll2008l) . The value 
we measure is consistent with the fractional luminosity/age re- 
lation from Decin et al. (2003). The disc extent, ~ 140 AU, is 
larger than other DUNES resolved discs, e.g. q' Eri (85 AU, 
ILiseau et al.ll2010l) or Ret (70-120 AU, Eiroa et al.ll2010h . but 
the derived black body temperature is also large, ~ 50 K, cf 60 K 
for the q' Eri disc which is around an F star and 30^0 K for the 
Ret disc which has a smaller extent around a star of similar 
spectral type. 

A complete analysis of the disc SED and physical structure, 
in which both classical power-law particle size distribution mod- 
els and a self-consistent collisional model are fitted to the avail- 
able photometric and imaging data will be presented by Lohne 
et al., 201 1 (in prep.). 



4. Conclusions 

We have presented Herschel PACS and SPIRE observations of 
HD 207129, the first to directly resolve the ring-like structure of 
a debris disc in both major and minor axes at far infrared wave- 
lengths. The disc extent (140 +32 AU), incHnation (51 + 5°) and 
typical (black body) temperature (~ 50 K), derived purely from 
the observations assuming no particular grain model, are similar 
to previous measurements. Compared to other discs around stars 
of similar spectral type and age, the disc of HD 207129 is both 
larger and warmer (for the observed size), though it is not com- 
pletely atypical in either respect, but situated at the margins of 
the known range of disc morphologies. A simple analysis based 
on a black body model has been proved unrealistic. A more de- 



tailed self-consistent study is left for future work (Lohne et al, in 
prep.). 
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Fig. 3. Original image (top), post stellar subtraction (middle) and Wiener deconvolved image (bottom) in the three PACS wave- 
lengths. The ring structure is clearly visible in both the 70 and 100 jim images. Image orientation is North up. East left. A distance 
scale bar has been provided, 10" is equivalent to 160 AU. The colour scale bar is in units of mJy/arcsec^ . The stellar position, 
derived from isophotes in the 70 and 100 observations, is marked as a white cross in each image. 
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